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ABSTRACT 

The Durham GALFORM semi-analytical galaxy formation model has been shown to re- 
produce the observed rest-frame 1500A luminosity function of galaxies well over the whole 
redshift range z = 5 — 10. We show that in this model, this galaxy population also emits 
enough ionizing photons to reionize the Universe by redshift z = 10, assuming a modest 
escape fraction of 20 per cent. The bulk of the ionizing photons is produced in faint galaxies 
during starbursts triggered by galaxy mergers. The bursts introduce a dispersion up to ~ 5 dex 
in galaxy ionizing luminosity at a given halo mass. Almost 90 per cent of the ionizing pho- 
tons emitted at z = 10 are from galaxies below the current observational detection limit at 
that redshift. Photo-ionization suppression of star formation in these galaxies is unlikely to 
affect this conclusion significantly, because the gas that fuels the starbursts has already cooled 
out of their host halos. The galaxies that dominate the ionizing emissivity at z = 10 are faint, 
with Mi5oo,ab ~ —16, have low star formation rates, M* ~ 0.06 /i _1 il/ yr _1 , and reside 
in halos of mass M ~ 10 9 /i _1 Af . 

Key words: galaxies: formation, galaxies: evolution, galaxies: high-redshift, intergalactic 
medium, cosmology: dark ages, reionization, first stars 
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1 INTRODUCTION 

Reionization refers to the transition in the state of the Universe from 
mostly neutral, following recombination at a redshift of z ~ 1000 , 
to highly i onized once more at later t imes. iGunn & Peterson! ill 96^) 
(and a lso iBahcall & Salpeterl 1 19651 ) realised as soon as I Schmidt! 
d 19651) published spectra of z ~ 2 quasars that the absence of 
significant Lyman-a absorption in their spectra implied that the 

2 ~ 2 Universe is very highly io nized. Tha t basic picture has 
not ch anged with the dis covery by iFan et al] d2003l) of z > 6 
QSOs ( iBecker et al]|2007t) or the novel method based on gamma- 
ray burst s as probes of the intergalactic medium (IGM) at e ven 
higher z l lTotani et al .11 20061 : IZafaret al]20ld : |Patel et all20ich . 

The fact that most of the hydrogen in the Universe is highly 
ionized at least as early as z ~ 7 is also consistent with the large 
Thomson scattering optical depth toward the surface of large scat- 
tering which is inferred from measurements of CMB fluctuations. 
This implies a 'reionization redshift' of 2 rc ion = 10.5 ± 1.2, if 
the transiti on from neutral to co mpletely ionized occurred instan- 
taneously dKomatsu et al j[2oToh - The temperature of the IGM de- 
pends on its reionization history because the thermal timescales 
are long: measurements of that temperature dSchave et alj|2000l) 
are also consistent with reionization occurring around z ~ 10 
dTheuns et a"Hl2002l) . 
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The current paradigm as to how reionization happens is that 
initially small HII regions form around individual sources of ion- 
izing photonsQ. As the sources become brighter and more nu- 
merous, isolated HII regions grow, merge, and eventually per- 
colate th roughout the IGM, see fo r example the early simula- 
tions by iGnedin & Ostrikeil dl997l) . The nature of the sources 
of the ionizing radiation is still unknown. While a number of 
works show that the majority of i onizing radiation i s proba- 
bly p r oduced by stellar sou r ces (e.g. Madau et a Jll999l:lGnedin] 
200C ; ISokasian et al.l 120031: ICiardi et all |2003l : iFurlanetto etal] 
2004 : Urac & CerJ 120071 : iTrac & Gnedinl l2009h . the exact con- 
trib ution of Population III stars or quasars is under debate 
(seelChoudhurv & Ferraral|2007l: IWvithe & Cenl|2007l ; lLoebll2009l : 
Volonteri & Gnedinll2009[|salvaterra et al.ll20ia for recent exam- 
ples). 

Depending on its spin temperature, the not-yet ionized HI dur- 
ing the epoch of reionization (EoR) could be detected in either 
emission or absorption in redshifted 21 -cm radiation, either in the 
form of a global step in the spectrum, or indeed probing the remain- 
ing neutral regions in a partly ioni zed IGM dMadau et al J 1 19971 : 
Ishaver et al ] | 19991 : iTozzi etakhood) . Because most plausible ion- 
izing sources will be highly clustered, the ionized bubbles could 



1 However, a strong background flux of higher-energy radiation, for ex am- 
ple X-rays from accreting black holes, may 'pre-reionize' the Universe iObl 



2 Raicevic et al. 



grow to be quite large, and the epoch where the IGM is 50% ionized 
may be best suited for direct detection with current and future ex- 
periments, such as LOFAF0, 21CM^B MWaQ, and eventually the 
The promise of a direct observational probe has stimulated 
considerable interest in the EoR, see recent reviews by for exam- 



plelBarkana & Loebl(l200l[) : ICiardi & Ferrard ( l2W5l) ; lLoebl ( l2006l) ; 
iTrac & Gnedinl ( 1200% . 

The recent installation of the Wide Field Camera 3 (WFC3)/IR 
on the Hubble Space Telescope has made it possible to search 
for 2 > 6 galaxies using the 'Lyman-break' drop-out technique, 
with a number of authors reporting the discovery of galaxies with 
z > 6 (based on their colours), with candidates up to z ~ 10 
teunker et al J 120091: IBouwens et alj|2007l. 12009 j ; iBouwens etall 
l2009d :l Oesch et al.l2009l). Are these the galaxi es that caused reion- 
ization? The analysis by iBunker et al.l d2009h suggests that these 
galaxies are unlikely to produce sufficient ionizing photons to 
reionize the Universe. In fact even at lower z ~ 6 there seems to 
be a problem, in the sense that the observed galaxies do not appear 
to produce sufficient phot ons to keep the IGM from recombining 
l lBolton & Haehneltll2007t) . 

In this paper w e use the GALFORM semi-analytic al model 
of galaxy formation dCole et alj|200b1;lBaugh et alj|2005l) to make 
theoretical predictions for the evolution of the emissivity of ion- 
izing photons from galaxies, e(z). The GALFORM model cal- 
culates the formation and evolution of galaxies in the framework 
of hierarchical structure formation in CDM, including baryonic 
physics such as gas cooling, star formation and supernova feed- 
back. In contrast to most other work modelling the contribution 
of galaxies to reionization, the GALFORM model which we use 
here was originally developed to try to explain the properties of 
galaxies at much lower redshifts. Predictions from GALFORM 
have been compared with a very wide range of observational data 
at lower redshifts z < 6, including galaxy luminosity functions, 
colours, sizes, morphologies, gas contents and metal li cities at red- 
shift z = dCole et al.l |200Ct iBaugh et alj 120051 ; iBower et al.l 
l2006t I Gonzalez et al.ll2009l) , the evolution of galaxies at optical , 
IR and sub-mm w a velengths feaugh et alj2005l ; lBower et al.l200r3 ; 
lLacev et al.ll2008l : iGonzalez-Perez et alj|2009h. and the evolution 
of Lj/a-emitting galaxies lEe Delliou et al J2006l : l0rsi et al.l2008h . 
In this paper, we use the lBaugh et al.l l l2005h varia nt of the GA L- 
FORM model. This model was already shown by IBaugh et alj to 
reproduce the observed rest-frame far-UV luminosity function of 
Lyman-break galaxies at z ~ 3, and the same model has recently 
been shown to reproduce the observed num bers of Lyman-brea k 
galaxies over the whole range z ~ 3-10 lEacev et alj|2010bb . 
Two important features of the iBaugh etalJ <2005l) model are that 
at high redshifts, most star formation happens in starbursts trig- 
gered by galaxy mergers, and the initial mass function (IMF) of the 
stars formed in such bursts is top heavy, containing a much larger 
proportion of high-mass stars than is found in more quiescent star 
formation environments such as the solar neig hbourhood. T his top- 
heavy IMF was introduced into the model bv lBaugh et alj in order 
to explain the observed numbers and redshifts of the faint sub-mm 
galaxies, now known to be very luminous, dust-obscured starbursts 
at z ~ 1 — 3. Models assuming a standard solar neighbourhood 
IMF were found to underpredict the numbers of sub-mm galaxies 



by an order of magnitude, if these models were also constrained to 
reproduce the present-day galaxy luminosity functions. 

In the present study, we want to investigate whether a model 
that is consistent with the new measurements of the Lyman-dropout 
galaxy population at z > 6 can produce sufficient ionizing photons 
to reionize the universe at z Te ion > 10. If so, we want to quantify 
which galaxies dominate e, and which aspects of the model affect e 
most. A si milar analysis based o n the GA LFORM model was p er- 
formed by iBenson et alj d2006h (see also iBenson et al.l l2002al lbh. 
who also included a simple model for the evolution of the HII 
volume filling factor. The present paper looks in more detail at 
the properties of the galaxies that cause reionization and the dark 
matter halos that host them, and how these are connected to the 
newly discovered z > 6 drop-outs. The properties of Lyman-break 
galaxies predicted by GALFORM over the whole redshift range 
z = 3 — 20 have been anal yzed in more detail in a companion 
paper by lLacev et alj fcOlObh . which also makes a more detailed 
comparison with the observed far-UV luminosity functions. We 
emphasize that the default values of the GALFORM model pa- 
ra meters used in the present work are identical to those chosen 
bv lBaughetafl lEool). which were adjusted to match a range of 
observed galaxy properties at lower redshifts. We will couple the 
GALFORM source model with a radiative transfer scheme to in- 
vestigate the progression of re ionization in more detail in a follow- 
up paper dRaicevic et alj|201oT) . 

This paper is organised as follows. In Section 2 we briefly 
discuss the main ingredients of GALFORM, paying particular at- 
tention to those aspects that most affect the ionizing luminosities of 
the galaxies. In Section 3 we show the evolution of the emissivity 
e for the default GALFORM parameters, discuss which galaxies 
dominate e, and how changes in GALFORM parameters affect e. 
In section 4 we show the corresponding far-UV luminosity func- 
tions at 2 = 6 and 10, explore the extent to which the currently 
detected galaxies constrain e, and how future observations with 
e.g. the James Webb Telescope will improve our understanding. We 
summarize in Section 5. 



2 METHOD 
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The GALFORM semi-analytical model fcole et alj kood) com- 
putes how galaxies form and evolve in the hierarchically growing 
dark matter halos of a cold dark matter Universe. The evolution of 
the halos themselves is described by halo merger trees, which are 
either extracted from an N-bo dy simulation or com puted using a 
Monte-Carlo scheme based on lLacev & Cole] i 19931) and improved 
bv lParkinson et alj j2008l) . The semi-analytical algorithm incorpo- 
rates physically motivated recipes for gas cooling, star formation, 
feedback from supernovae, galaxy mergers, metal enrichment, dust 
production and other processes, and in particular allows a calcula- 
tion of the observable properties of each galaxy, notably its broad- 
band l uminosity and colours, and its ionizing emissivity; see lBaugh] 
j2006|) for a recent review of semi-analytical methods. 

The buildup of dark matter halos of course depends on the 
assumed cosmological parameters, but the properties of the galax- 
ies associated with them are at least equally strongly depen- 
dent on the 'gastrophysics' governed by GALFORM parame- 
ters; for this reason we only consider th e cosmological param- 
eters used in the Millennium simulation ISpringel et alj d2005h . 
(fi m , Q A , fib, h, a 8 ,n s ) = (0.25, 0.75, 0.045, 0.73, 0.9, iff" 



Note that the oris inal lBaughetai]j2005l) model used a slightly different 
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Even at redshift z = onl y a very small fractio n of baryons 
have been converted into stars dFukugita et alj[l998l) . In particu- 
lar the faint- end slope of the z = K-band luminosity function, 
ql ~ — 1 dCole et all 1200 ll) . is much flatter than the low-mass 
slope of the dark halo mass function, qm ~ — 2 jPress & Schechterl 
[l974h . Therefore a crucial ingredient of any successful galaxy for- 
mation model is strong negative feedback to quench the fo rmation 
of small galaxies dWhite & Frenkll99ll ; lBe"nson et al.l2003h . GAL- 
FORM incorporates this and other effects with a set of rules, each 
with an associated set of parameters. Some of these have a large 
effect on the properties of early galaxies, others mostly affect the 
present-day galaxy population. Recent studies using GALFORM 
have c oncentrated on two different variants, that of Baugh et all 
d2005h (hereafter, BAUGH05) and of lBower et all J2006h (hereafter. 
BOWER06), which adopt somewhat differe nt prescriptions for star 
formation, feedback and the IMF (see also lLacev et all d2008l) for 
more details about the BAUGH 05 model). Th e BAU GH05 model 
includes superwinds (following iBenson et al. in order to 

better reproduce the bright-end of the optical and near-IR galaxy 
luminosity function at z = 0, while the BOWER06 model instead 
accompl ishes this by includin g feedback from accreting black holes 
(see also lCroton et al. I J2006I) ). The other most important difference 
between the two models is that the BAUGH05 model assumes that 
stars form with a top-heavy initial mass function (IMF) in star- 
bursts, and a normal solar neighbourhood IMF in galaxy discs, 
while the BOWER06 model instead assumes that all star forma- 
tion occurs with a solar neighbourhood IMF. In addition to this, the 
two models make somewhat different assumptions about the star 
formation timescale in discs, supernova feedback, the timescale for 
ejected gas to be re-incorporated into halos, and the triggering of 
starbursts. 

While the BAUGH05 and BOWER06 models predict similar 
galaxy luminosity functions at optical and near-IR wavelengths at 
z — 0, the BAUGH05 model is in much better agreement with the 
observed numbers of star-forming galaxies seen at high redshifts, 
selected either as Lyman-dropouts or f rom their sub-mm emission 
( iBaugh et alj|2005l : lLacev et alj|2010bl) . As we will show later, the 
BAUGH05 model also predicts higher ionizing emissivities at high 
redshifts than the BOWER06 model, and a correspondingly higher 
redshift of reionization, in better agreement with current observa- 
tional constraints. For these reasons, we concentrate in this paper 
on predictions from the BAUGH05 model. Not surprisingly, neither 
superwinds nor AGN greatly affect the predictions for galaxies at 
z > 6, since the massive galaxies that are affected by these pro- 
cesses are extremely rare at such early times. Nevertheless, we do 
find quite significant differences between these two popular GAL- 
FORM variants in what they predict for ionizing emissivities at 
z > 6, which are related to their different assumptions about star 
formation, supernova feedback and the IMF. We now discuss the 
physical processes incorporated in the BAUGH05 version of GAL- 
FORM model that have a large effect on z > 6 galaxies, and why 
they were introduced in the original model. 



2.1 Star formation 

The model assumes two distict modes of star formation, quies- 
cent star formation in galaxy discs, and starbursts triggered by 
galaxy mergers. In both cases the instantaneous star formation rate 
is parametrized as: 

Mcold 







(1) 



cosmology, (Q m , Q A , fi b , h, <r 8 , n s ) = (0.3, 0.7, 0.04, 0.7, 0.9, 1). The 
change of cosmological para meters was introduced for c onsistency with 
the Millennium-II simulation an-Kolchin et alj|2009l) which we will 
employ in future numerical simulations of reionization. 



where M co id is the amount of cold gas in the galaxy, and r* the star 
formation time scale. Neglecting the life-times of massive stars (the 
instantaneous recycling approximation), the stellar mass in long- 
lived stars then builds up at a rate 

M* = (1- J R)V ! (2) 

where R is the recycling fraction, see ICole et alj JioOOT) for more 
details. 

In the quiescent star formation mode, t* depends on the circu- 
lar velocity, Vdi sc , of the galactic disc at the half-mass radius, 
as t* = r*,o(Vdisc/200 kmr 1 )"*, with r*,o = 8 Gyr and 
a* = —3. This parametrization yields reasonable gas masses and 
star formation rates at low redshifts z ~ 0, and implies that ip is 
quite low at high redshifts. This makes the high-z discs gas rich, so 
that when galaxies me rge, there is a large reservoir of gas available 
for fueling a starburst l lBaugh et al. 12003) . 

Bursts of star formation are assumed to be triggered by galaxy 
mergers under certain conditions. The model includes both major 
and minor mergers, distinguished by the mass ratio of merging 
galaxies. Major mergers between spirals are assumed to destroy 
both discs and consume the remaining gas in a starburst. Minor 
m ergers were introduced in the model motivated by the simulations 
of Irlernquist & Mihosl dl995t) : such a merger does not destroy 
the disc, but does build up the bulge. The star formation time 
scale in the burst mode is shorter than in the quiescent mode (see 
IBaugh et alj|2005l) . 

The stellar initial mass function for quiescently forming stars 
is assumed to be similar to what is obse rved in the sola r 
neighbourhood, specifically that proposed by iKennicutl dl983l) . 
dN/dln(m) tx m~ x , with x = 0.4 for m < 1M and x = 1.5 
for m > lAf , However, in bursts the IMF is assumed to be 
top-heavy, x — 0. In both cases, the IMF covers the mass range 
0.15 < m/M Q < 120. 

Star formation with a top-heavy IMF in bursts triggered by 
gas-rich galaxy mergers results in large UV luminosities from the 
massive young stars, and also the production of large quantities of 
metals and dust from supernovae. This dust in turn absorbs the co- 
pious UV radiation and re-radiates it at far-IR wavelengths. Both 
the frequent bursts at high redshifts and the top-heavy IMF are 
needed to boost the number of very luminous high-z IR galax- 
ies to a level consistent with the observed number counts and 
redshift distribution of sub-mm galaxies. The parameters in the 
BAUGH05 model were chosen to match this sub-mm data, while 
at the same time yielding good fits to the Lyman-break galaxy lu- 
minosity function at z ~ 3, an d remaining conist ent with obser- 
vational constraints at z = dBaugh et alj 12003) . The case for 
a top-heavy IMF for the formation of at least a fraction of stars 
is further supported by the fact that its use during starbursts also 
results in better agreement with observed metallicities (including 
a/Fe ratios) in intracluster gas in clusters and stars in elliptical 
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galaxies dNagashima et al.l l2005tJlbh . Other independ ent observa- 
tional e vidence for variations in the IMF is discussed in lLacev et al.l 
We emphasize that our results do not depend crucially 
on the precise form of the top-heavy IMF assumed - similar re- 
sults would be obtained for an IMF in which the high-mass slope 
was fixed but the low mass turnover was varied, as proposed by 
iLarsorj Jl998l) . We will show below that the bursts, and the associ- 
ated change in the IMF during bursts, both have large effects on the 
emissivity of ionizing photons by GALFORM galaxies at z > 6. 



2.2 Supernova feedback 

The fact that galaxies in low-mass halos form stars very in- 
efficiently is likely due to energy injection from supemovae 
jPekel & Silkll 19861) . In the BAUGH05 model this is implemented 
by ejecting gas out of a galaxy disc at a rate 



Vdis, 



(3) 



so that it is no longer available for star formation. Here, Vdisc is the 
circular velocity of the galactic disc at the half -mass radius. Values 
of Vhot = 300 km s _1 and ahot = 2 were chosen to reproduce 
the fai nt-end slope of th e B-band galaxy luminosity function at 
z = teaugh et alj2005h . Such strong feedback also significantly 
quenches star formation in small halos at z > 6, and therefore 
has a large impact on reionization. Note that the BOWER06 model 
incorporates even stronger SN feedback in small halos. 

2.3 Photo-ionization feedback 

Star formation in small galaxies may be quenched as t he IGM be- 
come s ionized, either because cooling is suppressed dEfstafhioi] 
fl992h . or because the higher IG M gas pressure inhibits gas from 
falling into halos l lGnedin|[200ol) . or because photo-heating causes 
small galaxies to lose their gas l lHoeft et"al]|2006l : lOkamoto et al.l 
120081) . These effects may lead to a global suppression of star for- 
mation dur i ng and after the EOR, as seen in the simulations of 
ICrain et alj d2009l) . The standard approach in GALFORM is to 
model this by suppressing the cooling of halo gas onto the galaxy 
when the host halo circular velocity is below a threshold value 

Voire < Vcut ■ (4) 

at redshifts z < z cnt (but see also lBenson et alj d2002al) for a more 
detailed treatment). 

The default value of V cut = 60 k m s" 1 in the B AUGH05 
model, originally guided by the results of iGnedinl d2000h . is con- 
siderably larger than values found from more recent simulations 
dHoeft et alj|2006t lOkamoto et al]|2008h . The original BAUGH05 
model also assumed z cut — 6. Interestingly, because only the gas 
cooling in the halo is suppressed in GALFORM, a small galaxy 
with circular velocity V < Vcut can continue to form stars until 
it has exhausted its supply of cold (i.e. already cooled) gas. This 
way of suppressing galaxy formation in small halos once the IGM 
is ionized has consequences for reionization and also for the lumi- 
nosity function at later times, as we show below. 



2.4 Modeling Lyman-continuum and broad-band SEDs 

The GALFORM code computes the spectral energy distribution 
(SED) of each galaxy, given its star formation history and abun- 
dance evolution. The population synthesis models are based on the 



Padova stellar evolution track s combined with Kurucz model atmo- 
spheres dBressan et al.ll 19981). The dust extinction is modeled with 
a prescriptio n described by ICole et al. I d2000h with improvements 
described in lLacev et alj fcOlObl) . Convolving the SED with a fil- 
ter response yields broad-band luminosities for the galaxy. Below 
we will use the rest-frame 1500A broad-band AB magnitudes of 
GALFORM galaxies to compare against observed galaxy lumi- 
nosity functions at approximately the same rest-frame wavelength, 
after reseating observed luminosities and number densities to the 
same MILLENNIUM cosmology as assumed in the model. 

GALFORM also computes the Lyman-continuum luminosity 
for each galaxy, expressed as the emission rate of ionizing photons, 

A^LyC — / t~~ dv , (5) 

J "thresh 

where L v is the SED of the galaxy and fthrcsh is the Lyman-limit 
frequency, ftfthreah = 13-6 eV. Note that the number of ioniz- 
ing photons produced per solar mass of stars formed is very differ- 
ent for the Kennicutt IMF assumed during quiescent star formation 
compared to the top-heavy IMF in bursts (A/ Ly c /Af* = 3.2 x 10 60 
and 3.5 x 10 respectively, for solar metallicity). 

A considerable fraction of those ionizing photons may be 
absorbed locally in the interstellar medium of the galaxy or by 
gas in the surrounding halo, and the fraction / osc of photons that 
does manage to escape into the IGM is very uncert ain. Observa- 
tions o f z ry 3 — 4 Lyman-brea k galaxies (LBGs) bv lSteidel et al.l 
j200lh andlShaplev et al.l J2006T> suggest / csc ~ 0.01 - 0.1 or even 
lower dGiallongo et alj|2002h (but note the slightly different def- 
inition of / osc there). The escape fraction may depend strongly 
on the porosity of the interstellar medium wi thin the galaxy or 
the presence of sup ernova-driven winds (e.g. ICiardi et aU 120021 : 
IClarke & Ovfosxm . Some of the more recent models that attempt 
to include these effects suggest that z > 6 galaxies may have sig- 
nificantly larger escape fractions, / esc ~ 0.5, (e.g. IWise & Cenl 
120091 : iRazoumov & Sommer-Lars en 2010). Note that this param- 
eter is unlikely to be independent of metallicity, gas content, and 
halo mass. In this paper, we simply assume / csc to be the same for 
all galaxies. 

We will now discuss the net emissivity of ionizing photons 
in the BAUGH05 model, and how that depends on GALFORM 
parameters. 



3 IONIZING EMISSIVITIES 

The emissivity e(z), the number of ionizing photons produced per 
unit comoving volume at redshift z, is found by summing the 
Lyman-continuum luminosity of all galaxies, per unit volume, 



(*) = 



iV LyC -I'(iV Ly c)diV' Ly c, 



(6) 



where <&(Nhyc) is the Lyman-continuum luminosity function. The 
emissivity e(z) increases by approximately 1.5 dex between z — 
13 and z = 5 in the BAUGH05 model (FigQ] thin line), mostly as 
a consequence of evolution in the halo mass function, as we will 
show below. 

Integrating e(z) down to a given redshift yields the total num- 
ber of ionizing photons produced per unit comoving volume up to 
that time. This number can be compared to the mean comoving 
number density of hydrogen atoms, tih. Reionization will occur 
when their ratio 
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Figure 1. The ratio TZ(z) of the number of ionizing photons produced per 
hydrogen atom up to redshift z in the two fiducial GALFORM models, 
BAUGH05 and BOWER06 (thick lines, left y-axis) as well as the total emis- 
sivity, e(z), in the same models (thin lines, right y-axis). The horizontal 
dashed lines mark the minimum number of photons per H atom that must 
be produced to achieve reionization: in the most optimistic case, only one 
(bottom line), but 10 or more when reasonable values for the ionizing es- 
cape fraction and mean number of recombinations per H atom are taken 
into account (top line). The BAUGH05 model produces ~ 100 times more 
ionizing photons at z ~ 10 than BOWER06 and reaches 10 photons per H 
atom Az ~ 5 earlier. The decreased slope in e(z) at z ^ 6 is caused by the 
turn-on of photo-ionization feedback at z = 6 in both models. 




Figure 2. Lyman-continuum photon luminosity, N^ydM, z), of halos as 
a function of halo mass M, in the BAUGH05 model at z = 10 (me- 
dian and mean relation are shown as thick and thin solid lines, respec- 
tively). N-^yQ increases approximately as A^yC oc M 1,8 for small halos 
M < 2 X 10 9 /i _1 M0, and as N^yC tx M for more massive halos, with 
little dependence on redshift. The 50 and 90 per cent ranges of N-^ y c(M) 
at given halo mass are shaded red and purple, respectively. There is up to 
5 dex range in Nl y c at a given mass, a consequence of the dominance of 
starbursts in producing ionizing photons. 



, , (" e(z) dz 
TZ(z) = Jo ° v ; — , (7) 
n H 

is TZ = (1 + iV r ec)//esc. Here, N IC c denotes the mean number of 
recombinations per hydrogen atom up to reionization, and / esc is 
the mean escape fraction from Section [2l4l 

Estimating N rcc is not straightforward. Recombinations can 
occur in the higher-density regions of the general IGM, in 'mini- 
halos' that have too shallow potentia l wells for star formation 
dShapiro et al]|2004l : ICiardi et al.|[200r3) . or in even higher-density 
regions associated with Lyman-limit or damped Lyman-a systems. 
The value of N rcc will itself depend on e(z) dz, since a slower 
build-up of the ionization rate will allow more time for recombi- 
nations. Interestingly, once the IGM is ionized, the smoothing of 
the density field due to gas pressure following photo-heating re- 
duces the recombination rate <Pawliketal]|2009l) . Current simula- 
tions of the EoR suggest values of N re c of a few jlliev et al]|2006l : 
iMcQuinn et al.ll2007UTrac & Cenll2007T) . 

Combining the estimate of 1 + N rcc ~ 2 with a reasonable es- 
cape fraction of / esc ~ 0.2 then suggests that reionization requires 
a value of TZ ~ 10. This is plotted as a function of redshift for 
the default values of the BAUGH05 and BOWER06 GALFORM 
parameters in Fig. [T] (thick lines), suggesting that the BAUGH05 
model will produce a reasonable reionization redshift z rc i on ~ 10, 
Az ~ 5 before BOWER06. Next we discuss the properties of the 
galaxies and halos that dominate the emissivity in the BAUGH05 
model, and how strongly these depend on the assumed parametriza- 
tion in the model, following the same order as in the previous Sec- 
tion|2] 




Figure 3. Main panel: Lyman-continuum emissivity as a function of halo 
mass, de(M, z)/d log 10 (M), for various redshifts indicated in the panel. 
The emissivity is low for very low-mass halos that are unable to cool gas, 
reaches a peak which increases with decreasing z, and a tail towards larger 
masses set by the exponential drop in the number of massive halos. At 
z ~ 10 most ionizing photons are produced by halos in a relative small 
mass range, ~ 1 dex. Top inset: cumulative fraction f c of ionizing pho- 
tons produced in halos more massive or less massive than a given value 
(rising and falling curves, respectively). The mass of halos below which 50 
per cent of ionizing photons is produced rises by approximately an order of 
magnitude from ~ 8 X 10 8 h^ 1 M Q at z = 14 to ~ 8 X 10 9 /i _1 M at 
z = 6. 
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Figure 4. Dependence of the total number of ionizing photons produced 
per hydrogen atom up to redshift z, IZ(z), on the starburst parameters in 
BAUGH05: default model (black), no bursts (red), including bursts, but not 
the change to a top-heavy IMF in bursts (blue). Including bursts increases 
e(z) by a factor 5-10, depending on redshift. The effect of the change in 
IMF in the bursts is large, yet even without it bursts still increase e by a 
factor of ~ 2. Neglecting bursts delays reionization (TZ = 10) by Az ~ 4. 
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Figure 5. Dependence of emissivity as a function of halo mass, 
de/dlog 1( ) M, on the burst parameters in the BAUGH05 model. The char- 
acteristic halo mass at which 50 per cent of the ionizing photons is pro- 
duced does not greatly depend on the burst parameters. However, switching 
off the bursts (red short dashed line) extends the halo mass range in which 
the majority (~ 90 per cent) of ionizing photons is produced by ~ 1 order 
magnitude in comparison to the default model (solid black line). 



3.1 Effect of star formation parameters and IMF 

The number of ionizing photons produced per unit time by galaxies 
in a halo of given mass, Al y c (M, z), is plotted as a function of M 
in Fig. [2] The virial temperature T V - 1T of halos with M < M m - ln ~ 
10 8 /i _1 Mq is too low to enable radiative cooling by atomic lines 
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Figure 6. Distribution of Lyman-continuum photon luminosities, A^LyC at 
z = 10, for halos with mass M fa 10 9 /i — x Mq. Different line styles refer 
to different models for the bursts, vertical dotted and solid lines indicate me- 
dian and mean Nh y c m the default model, respectively. The distribution of 
NhyC peaks at a few times 10 50 h photons s — 1 , but allowing bursts in- 
troduces a long tail towards much more luminous galaxies (red versus black 
histograms), with the change in IMF in bursts having a large contribution to 
this (blue versus black histograms). This tail makes the mean A^L y c nearly 
2 dex brighter than the median. In the default model with a top-heavy IMF 
in bursts there is a nearly 5 dex range in luminosity at given halo mass. 

and hence such halos do not form starfl Given that T v - lr oc (1 + 
z) at fixed M, there is strong redshift dependence in N^ y c{M, z) 
at very low masses, but above this minimum mass GALFORM 
predicts essentially no evolution in the mean A^LyC (M, z) between 
2 = 15 and 2 = 6, but with a modest ~ 50% decrease in the 
median in halos with mass M > 10 10 /i _1 Mq in the same redshift 
range. 

The mean A^LyC at a given halo mass increases approximately 
as iV Ly c oc M 1 - 8 for small halos M < 2 X 10 9 /i _1 M , and 
roughly as A^yC oc M for more massive halos, in contrast to many 
recent simulations of reio nization which assume a simple NhyC 
M relation for all M (e.g. lFurlanetto et alj2004llliev et alfeOOfj) . 
Interestingly, there is a very large difference between the mean and 
median of Ni, y c at given M, and there is also a very large range, 
up to ~ 5 dex, in A^yC at given M (Fig.|2]l- Both are consequences 
of the importance of bursts in generating ionizing photons, as we 
will discuss in more detail below. 

The total Lyman-continuum emissivity per dex in halo mass 
de/dlog 10 (Af ) (Fig.[3j, can be obtained by combining the mean 
luminosity of a single halo of given mass, ALyc(Af), with the 
number of halos of that mass, dn/dlog 10 (M). This function 
evolves rapidly as a consequence of the rapid build-up of more 
massive halos as time progresses. The halo mass below which 
50 per cent of ionizing photons are produced increases from ~ 
8 x 10 8 /i _1 M Q at 2 = 14 by an order of magnitude to ~ 
8 x 10 9 h' 1 M e at 2 = 6 (top panel of Fig. HJ. At high 2, the 



7 We recall that this GALFORM model does not consider Pop. Ill stars 
that form due to molecular cooling in such small halos. 
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mass range of halos that contribute significantly to e is relatively 
small, of order 1 dex, since it is limited at low M by M m i n and at 
large M by the exponential drop in the abundance of more massive 
halos. At later redshift z ~ 6, de/dlogioM is nearly independent 
of M over nearly 2 dex, a consequence of the fact that the ionizing 
photon luminosity of halos increases with halo mass approximately 
as N-Lyc(M) oc M 1 (dotted line in Fig. [2}, whereas the number 
density of halos decreases with increasing mass approximately as 
dn/dlog 10 M oc A/ -1 . 

The impact of starbursts on the emissivity is quantified in 
Fig. [4] In the default BAUGH05 model, bursts increase the ion- 
izing emissivity relative to that from quiescent galaxies both as a 
consequence of the reduction in star formation timescale, Eq. Q}, 
and because of the assumed change to a top-heavy IMF. The net 
effect is a factor 5-10 increase in e depending on redshift, with ap- 
proximately 65 per cent of the increase due to bursts following a 
minor merger. Most of the increase in A^LyC is a consequence of 
the assumed change in IMF. 

Neglecting bursts does not affect the 'characteristic' halo mass 
below which 50 per cent of the ionizing photons are produced 
(Fig. O but it does increase the range of halo masses responsible 
for the majority (e.g. 90 per cent) of ionizing photon production by 
~ 1 dex (compare solid black and short dashed red lines in the top 
inset of the same panel). 

Bursts skew the distribution of Nt, y c at given halo mass by 
introducing a long tail of much more luminous galaxies which hap- 
pen to be bursting, with again the assumed change in IMF playing 
a dominant role (Fig. [6j. These few, but relatively bright, galaxies 
dominate the emissivity at that halo mass by a large factor. Re- 
markably, there can be nearly a 5 dex range in Lyman-continuum 
luminosity at given halo mass. 

We conclude that bursts are a crucial ingredient in order for 
the BAUGH05 model to produce that many ionizing photons by 
2 ~ 10. Not only do stars form at a greater rate due to the de- 
crease in the star formation timescale, but especially the change to 
a top-heavy IMF in bursts, originally introduced to produce suf- 
ficiently luminous sub-mm galaxies at z — 1-3, and to produce 
sufficient metals by z = 0, causes a small fraction of galaxies to 
emit copious ionizing radiation. The bursts occur mostly due to mi- 
nor mergers, and are so effective because the merging galaxies are 
very gas rich, itself a consequence of the inefficient star formation 
in their quiescent state. Bursts also introduce nearly 5 dex of scatter 
in the A^LyC-halo mass relation. These same bursts are also a cru- 
cial ingredient for reproducing the observed l uminosity function of 
Lyman-break galaxies at z > 6, as shown in lLacev et al.l ( l2010bl) 
and also discussed below (Fig.[5). But first we investigate the effect 
of the feedback parameters on e. 



3.2 Effect of supernova feedback parameters 

We consider two variants to the default BAUGH05 supernova feed- 
back parametrization to investigate how strongly they affect the 
emissivity of ionizing photons. The 'weak' feedback choice, shown 
in Fig. UJ (green dashed line), uses parameters (Vhot> Qfoot) = 
(100 km s , 1) (as defined in Eq. 0, as opposed to the default 
BAUGH05 values of (300 km s , 2). The ionizing emissivity of 
the weak feedback model is not very different from a model without 
any SN feedback at all; it produces nearly twice as many ionizing 
photons as the default BAUGH05 model, increasing the reionization 
redshift, for which 1Z = 10, by Az ~ 0.7. The 'strong' feedback 
model has (Vhot,Qhot) = (500km s ,3), close to the values 
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Figure 7. Dependence of the total number of ionizing photons produced 
per hydrogen atom up to redshift z, 1Z(z), on the GALFORM parameters 
that govern supernova feedback. Decreasing the efficiency of SN feedback 
(green dashed line) doubles the production of ionizing photons, resulting 
in reionization occurring Az ~ 0.7 earlier. Increasing feedback from su- 
pernovae to the values used in the BOWER06 model (red dot-dashed line) 
delays reionization by Az ~ 2. 

(475 km s" 1 , 3.2) used in BOWER06 ; this choice of parameters 
decreases e(z) by a factor ~ 5, delaying reionization by Az ~ 2. 

Even stronger feedback is probably ruled out by the compari- 
son with the observed z — 6 Lyman-break far-UV LF discussed in 
Fig.[T0]below, but all three models are probably equally consistent 
with the z = 10 LF. This is not surprising since the SN parameters 
affect mostly the fainter galaxies that are currently below the de- 
tection limits at these very high redshifts. We note that the standard 
approach in GALFORM modelling is to constrain the SN feed- 
back parameters by comparison with galaxy properties at z = 0. 
However, even if one chooses to relax the z = constraints on 
the SN feedback, on the grounds that SN feedback might operate 
differently in early galaxies, the constraints on this from the 2^6 
Lyman-break LFs still limit the uncertainty in e to a factor ~ 2 in 
the BAUGH05 model. 

3.3 Effect of photo-ionization feedback parameters 

As discussed in Section [3~3l the effect of photo-ionization feedback 
from reionization on galaxy formation is modeled in GALFORM 
with a simple prescription, whereby gas cooling is suppressed in 
all halos of circular velocity Kirc < V cu t after the reionization 
redshift z cut , Eq. {4}. The key feature of this prescription is that 
the cold gas already present in galaxies before the onset of photo- 
ionization feedback is allowed to form stars after z cut . This results 
in a significant delay between the time at which the surroundings 
of the galaxy become ionized and the quenching of star formation. 
This is in contrast to several current simulations of re ionization, 
which assume that suppression is instantaneous (e.g. Illiev et all 
l2006t) . The delay is in fact so large that the suppression of star for- 
mation (and hence also the production of ionizing photons) due to 
photo-ionization has little effect on the progression of reionization, 
as we will show elsewhere. 

However, given enough time, photo-ionizing feedback does 
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Figure 8. Dependence of Lyman-continuum emissivity as a function of halo mass, de/dlog 10 M, at redshift z = 5, on the GALFORM parameters 
that describe photo-ionization feedback. Left panel: dependence on the reionization redshift z cu t (indicated in the legend; 'no z cxl t assumes reionization 
occurs below z = 5) below which the IGM is assumed to be fully ionized. The impact of photo-ionization suppression takes a long time to take effect, 
but when suppression sets in it dramatically reduces the LC luminosity of the galaxies. Right panel: dependence on the circular velocity V cu t below which 
photoionization feedback affects the galaxy. The suppression in de/dlog 10 M becomes evident at circular velocities below ~ 2V cu t. The scale on top of 
both panels gives the circular velocity of the halos (in kms -1 ) at z = 5. The numbers next to the lines give the ratio of the total emissivity of that model 
compared to the model with no photoionization feedback (black lines). 



have a strong effect on the ionizing emissivity, as shown in Fig. [8] 
Note that the default BAUGH05 model uses a value of V cu t = 
60 km s _1 which is unrealistically high compared to more re- 
cent simulation results, which reduces e by as much as 50 per 
cent by redshift 5 compared to the no reionization model (assum- 
ing reionization occurs at z cut — 10). The more modern value of 
V cu t ~ 30 km s _1 , suggested by the simulations of lOkamoto et al.l 
yields a smaller yet still significant decrease in the total 
emissivity at z = 5 of 15 per cent. 

We conclude that photo-ionization suppression as imple- 
mented in GALFORM has little effect on the production of ioniz- 
ing photons until well after reionization, but it does affect the emis- 
sivity at later times. Interestingly, the photo-ionization quenching 
of star formation also has observable effects on the Lyman-break 
LF, as we discuss in more detail below (Fig. II It. 



4 FAR-UV LUMINOSITY FUNCTIONS OF THE 
GALAXIES THAT CAUSED REIONIZATION 

The Lyman-break colour-selection technique has proven to be very 
effective for identifying large samples of star-formin g galaxies at 
high redshifts since its first application at z ~ 3 dSteidel et al.l 
19961). This selection method was first applied at 2 ~ 6 by 
Bouwens et alj d2003T) . and recent deep near-IR imaging with Hub- 
ble Space Telescope (HST) has been used to discover significant 
numbers of candidate Lyman-b reak galaxies (LBGs) at z ~ 7 — 8, 
and a few candidates at 2 ~ 1 llBouwens et alj2007l:lBunker et al 



20091 : iBouwens et all l2009d : iBouwens et alj l2009d : Toesch et al 



2009). We therefore now have direct detections of a part of the 
galaxy population responsible for r eionizing the univers e at z ~ 
6 — 10. The companion paper by lLacev et"aL1 ( feOlObl) presents 



a detailed comparison of the predictions of GALFORM models 
with observations of Lyman-break galaxies over the whole redshift 
range 2 = 3 — 10, including rest-frame far-UV luminosity func- 
tions, sizes, masses and other properties. In this section, we inves- 
tigate what constraints can be put on the GALFORM parameters 
to which the emissivity of ionizing photons e is particularly sensi- 
tive from observations of the rest-frame far-UV (1500A) luminosity 
functions of 2 ~ 6 — 10 Lyman-break galaxies alone. We also in- 
vestigate the extent to which the currently observed Lyman-break 
galaxies contribute to the total emissivity of ionizing photons, ac- 
cording to the GALFORM model. 



4.1 Effect of star formation parameters and IMF 

The rest-frame 1500A broad-band GALFORM LFs at 2 = 6 and 
2 = 10 are compared against the HST data on LBGs in Fig. [9] The 
default BAUGH05 model reproduces the LFs at both redshifts, a 
considerable success. Clearly, starbursts are crucial for bringing the 
1500 A luminosities of the galaxies to the observed levels (compare 
the red short dashed lines for the model without bursts with the 
other two lines). These same bursts also produce the bulk of the 
ionizing photons, as we showed in Fig. [5] 

Interestingly both the model with a top-heavy IMF in bursts 
(the defau l t mod el, black lines), and a model which uses the same 
iKennicuttl Jl983l) IMF in both quiescent galaxies and bursts (blue 
dashed lines) fit the observed LFs nearly equally well at these red- 
shifts, notwithstanding the significant differences between these 
models that we pointed out in, for example, Fig. [4] The reason for 
this is dust extinction: the default model with the top heavy IMF 
produces more m etals and hence also more dust as compared to the 
IKennicuttl Jl983h IMF, and the larger dust extinction partly com- 
pensates the larger intrinsic far-UV luminosities (see lLacev et alj 
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Figure 9. Rest-frame 1500 A broad-band luminosi ty functions of t he de - 
faul t BAUGH05 model (line s) compared to data from lBouwens et al] j2007l) 
and iBouwens et af] fe009al) . at redshifts z = 6 and 10 (symbols with er- 
ror bars; downward pointing arrows mark Icr upper limits). Both the de- 
fault BAUGH05 model (black solid lines) and the single IMF variant (long 
dashed blue lines) produce reasonable fits to the observed LFs at both red- 
shifts. The insets in each panel show the cumulative fraction of ionizing 
photons produced in galaxies brighter than, or fainter than, a given value of 
the Afi5oo,AB absolute AB magnitude (rising and falling curves, respec- 
tively). 

2010b. for more details). Previously we found that a change in IMF 
affected the ionizing emissivity considerably (Fig. [4}, but there we 
assumed that the escape fraction of ionizing photons / csc , is sim- 
ply a constant. A physically motivated / csc would presumably de- 
pend on galactic dust content, reducing the differenc e between the 
top-h eavy IMF and single IMF emissivities (see e.g. iBenson et al .1 
which would shift the completion of reionization we found 



Figure 10. The effect of the supernova feedback parameters on the pre- 
dicted rest-frame 1500 A luminosity functions in the BAUGH05 model at 
redshifts 6 (top) and 10 (bottom), and models with weaker and stronger 
feedback (green and red lines, respectively); the corresponding emissivities 
were shown Fig. [7] The data (solid points) are from Bouwens et al., as in 
Fig|9] The weak feedback model (green dashed line) slightly over predicts 
the number of galaxies at z ~ 6, and the strong feedback model under pre- 
dicts the numbers. However at z ~ 10 the bright, observed end of the LF is 
equally well fit by all models. 

here to lower redshifts. We will examine these issues in future 
work. 

The currently detected candidate LBGs contribute only a 
small fraction of the total emissivity of the whole population of 
galaxies predicted by GALFORM at high-z. Even at z ~ 6 
(top panel), galaxies brighter than the current observational limit 
(Afi5oo,AB,min ~ —18) contribute only ~ 40 per cent of the to- 
tal ionizing emissivity (solid black line in the top inset). If a single 



10 Raicevic et al. 



m 1 

< 



VA 





p. 



bfi 

-a 



o 



-3 
-4 
-5 

-6 















■ * 1 i' 1 I i i i i i i i i" - "" i — 


"I T 




i i i i i i i i i i 






* ML* ™ ''/, X 


z=6 : 






• Bouwens obs. \ 






; default (z cut =6) 






" z cut =1 5 V cut= 6 W s 






! z cut= 12 > V cut= 60 W s 






: z cut= 10 > V cut= 30 W S 




i . ; 



-14 



-16 



-18 

M 



-20 



-22 



1500 AB 



Figure 11. The effect of photo-ionization on the predicted rest-frame 
1500 A LF in the BAUGH05 model at redshift 6. The models differ in their 
choice of reionization redshifts (z cut ), and of the halo circular velocity be- 
low which galaxies are affected by photo-ionizing feedback (Vcut). The 
corresponding emissivities were shown in Fig[8] The data (solid points) are 
from Bouwens et at, as in Fig[9] If galaxies with V C ut=60 km s _1 are af- 
fected by suppression, then early reionization (z C ut ^ 10) can be ruled out 
by the current data, since then the predicted number density of galaxies at 
-^1500, AB ~ —18 is ~ 4 times lower than observed (red lines). A more 
reasonable suppression scale of V^ u t=30 km s _1 is consistent with early 
reionization (green lines). 



lKennicuttlfl983T) IMF is assumed, that fraction is even lower (~ 20 
per cent; long dashed blue line). At z ~ 10 (bottom panel), more 
than 90 per cent of ionizing photons are emitted by galaxies below 
the current detection limit for the default BAUGH05 parameters, 
and for a single IMF model that fraction is ~ 95 per cent. 

The BAUGH05 model predicts that the galaxies that produce 
the bulk of the ionizing photons at z ~ 10 are intrinsically faint, 
with 50 per cent of ionizing photons produced in galaxies fainter 
than 77iab ~ 31 in the H-band. Clearly it will be challenging to de- 
tect a significant fraction of the galaxies that emit the photons that 
reionized the Universe, even with the James Webb Space Telescope, 
see e.g. the JWST white paper by Stiavelli et a/Jfl 

4.2 Effect of supernova feedback parameters 

The strength of supernova feedback cannot be strongly constraine d 
with the current z > 6 data (Fig|T0j see also lLacev et alj|2010bh . 
At the lowest redshift (z = 6; top panel), the faint end currently 
probed provides some constraints on the strength of the supernova 
feedback, with the weak and strong models on either side of the 
data. However, the z = 10 data only probes the very brightest 
galaxies, for which all three models predict very similar LFs. 

Of course, the supernova feedback parameters in GALFORM 
are strongly constrained by even lower redshift data. However, the 
reader should keep in mind that the emissivities we predict here 

8 http://www.stsci.edu/jwst/science/whitepapers/ 



are contingent on the assumption that the basic physics of galaxy 
formation (in particular the impact of supernova feedback on regu- 
lating star formation) is the same at all redshifts. If for some reason 
this is not true, the currently available observations at z > 6 do not 
probe sufficiently faint galaxies to determine the impact of super- 
nova feedback on the total emissivity produced by all galaxies. 



4.3 Effect of photo-ionization feedback parameters 

As discussed in Section [331 the high-z 1500 A LF may hold in- 
formation about the reionization history, if star formation in galax- 
ies is quenched once their surroundings are ionized. This is illus- 
trated in Fig. [TT] The z — 6 LF is reasonably well fit by the de- 
fault BAUGH05 model, which assumes that reionization occurs at 
^cut = 6 (and hence for which there is no suppression in Fig. II It. 

However, recent CMB measurements of the Thomson scatter- 
ing optical depth suggest reioni zation at z ~ 10, ass uming an in- 
stantaneous reionization model iKomatsu et al]|201fj|) . The GAL- 
FORM model with such early reionization and Vcut = 60 km s - 
underpredicts the faint end of the observed z = 6 luminosity func- 
tion by a considerable amount, a factor ~ 4 for galaxies with 
A^i500,ab fainter than -18. Clearly photo-ionization suppression is 
then too strong. But we already argued that the default value of the 
halo circular velocity below which galaxies are affected by photo- 
ionizing feedback (V cu t = 60 km s -1 ) is too hi gh, with the hydro- 
dynamical simulations o f lOkamoto etai] J2008h suggesting a much 
lower value of V cu t = 30 km s _ . With this lower value of Vcut, 
the LF at z = 6 is in good agreement with the data, even for an 
early re ionization redshift (green dotted line, see also lLacev et all 
in fact this model fits the z = 6 data best. Noting that the 
CMB data is the strongest current constraint on reionization, we 
argue that this result gives an observational constraint on the char- 
acteristic strength of photo-ionization feedback that strengthens the 
conclusion from current simulations. 

The far-UV luminosity fun ctions predicted by th e BAUGH05 
model and presented here and in lLacev et alj J2010bh show a very 
good a greement with the z > 6 data o flBouwens etalj ( l200ll200a . 
12009 al ibi). This is a significant success for a model for which the 
parameters were chosen to match much lower redshift data, and 
provides us with reasonable confidence in using the ionizing lumi- 
nosities predicted by this model in future, more detailed modeling 
of the reionization process l lRaicevic et alj201ol) . 

We have seen that the BAUGH05 model predicts that the bulk 
of ionizing photons is produced by galaxies significantly below the 
current detection limit. It is a common practice to fit observed LFs 
with a Schechter function, and use the fit to extrapolate the LF to 
fainter galaxies. We show in the Appendix that this approach can 
lead to significant errors in estimating the total emissivity, since the 
LFs predicted by GALFORM deviate significantly from Schechter 
functions in some ranges of luminosity, in particular due to the ef- 
fects of bursts. As a result, the Schechter fit parameters depend sig- 
nificantly on the luminosity range over which the fit is done, and 
the total emissivity estimated by extrapolating this fit is sensitive to 
the minimum luminosity set by the observational detection limit. 



5 CONCLUSIONS 

We used the iBaugh et alj J2005h version of the GALFORM 
galaxy formation model to compute the emissivity (e) of hydrogen- 
ionizing photons in the redshift range relevant for reionization, 
z > 6, and investigated the impact of changing some of the model 
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parameters from their default values. A crucial element of this 
model is that mergers between gas-rich galaxies increase e dra- 
matically compared to a model without bursts, mainly due to the 
change to a top-heavy IMF in bursts assumed in the model. The 
iBaugh et al.l model, with the same parameter values as used here, 
has previously been shown to reproduce a wide range of observed 
galaxy properties at lower redshifts. 

The main points presented in the paper are: 

• The BAUGH05 model produces enough ionizing photons to 
complete reionization by 2 ~ 10 with galaxies alone, assuming 
a reasonable photon consumption (2 photons per hydrogen atom, 
allowing an average of 1 recombination per H atom) and a 20 per 
cent escape fraction of LC photons from galaxies (Fig.Q~|l. 

• Starbursts are crucial for boosting the ionizing emissivity lead- 
ing up to reionization. The majority of ionizing photons is produced 
in a relatively small fraction of galaxies at any given time that are 
bursting, and that are up to 5 dex brighter than non-bursting galax- 
ies in halos of the same mass. Such bursts also increase the impor- 
tance of intermediate-mass halos (M ~ 1O 9 /i -1 M0) compared to 
simpler models that do not include bursts (Fig.[4]l. 

• The top-heavy IMF used in the burst star formation mode is 
the main factor making the bursts so luminous, with ~ 10 times as 
many ionizing photons emitted p er solar mass of stars formed as 
compared to the lKennicuttl (l983) IMF. The change to a top-heavy 
IMF in starbursts was previously introduced in the model to repro- 
duce the sub-mm galaxy counts at lower redshifts (z ~ 1 — 3), not 
the ionizing emissivity we discuss here, but it is crucial for com- 
pleting reionization in agreement with current observational con- 
straints. The model with a single IMF reionizes Az ~ 2.5 later 
than the default model (Fig.|4). 

• The assumed strength of supernova feedback has a strong im- 
pact on the ionizing emissivity, because the galaxies that dominate 
e reside in relatively low-mass halos (Fig.|7J. This fact is of course 
well known at lower redshifts where strong feedback is required 
to reproduce the faint-end of the galaxy luminosity function (e.g. 
IColeetalJl200Ch . but is often ignored in reionization modeling, 
where a simple linear mass-luminosity relation is assu med. 

• As also shown in the companion paper by lLacev et"al] 
l l2010bh . the BAUGH05 model reproduces the observed z ~ 6 — 10 
rest-frame 1500A luminosity functions well (Fig. |9), with bursts 
a crucial ingredient in boosting the UV luminosities of galaxies 
to the observed levels. The good agreement between the predicted 
and observed UV luminosity functions gives credence to using the 
model for computing e as well. In the model, ~ 90 per cent of ion- 
izing photons are produced by galaxies that are below the current 
HST detection limit at z = 10, with 50 per cent of ionizing photons 
produced by galaxies fainter than ibab ~ 31 in the H-band. The 
intrinsic faintness of the sources will make it very challenging to 
detect a significant fraction of the galaxies that caused reionization, 
even with JWST. 

• The shape of the rest-frame far-UV luminosity function in the 
BAUGH05 model resembles a Schechter function, but with signifi- 
cant departures due to bursts. Given that the z > 6 data only probe 
the bright end of this LF, extrapolating a Schechter function fit to 
estimate the contribution from galaxies below the detection limit 
can be inaccurate (see the Appendix). 

As in all models of reionization, a significant uncertainty is 
the fraction / csc of ionizing photons produced by galaxies that can 
actually escape into the IGM. We have intentionally used a simple 
estimate for / csc , and our default value of 20 per cent is somewhat 
higher than found observationally in lower redshift observational 



studie s (e.g. / osc ~ 10% for LBGs at z = 3 - 4 ISteidel et all 
l200ll) . A high dust content, one of the consequences of using a 
top-heavy IMF, may d ecrease the escape fraction by as much as 
an order of magnitude l lBenson et all2006h . On the other hand, the 
fraction of the ionizing photons that can escape into the IGM dur- 
ing a burst could be significantly increased over the escape fraction 
during quiescent star formation, due to the galactic wind driven by 
the starburst. Detailed numerical models that include turbulent mo- 
tions of gas in small galaxies find that f esc can be as high as 0.5 - 1 
during a burst jFuiita et alj|2003l: Iwise & Abelll20o3 ; Iwise & CerJ 
120091 ; iRazoumov & Sommer-Larsenl 2010l) . The enhancement of 
/esc in bursts is likely to be more dramatic for smaller galaxies 
than for larger ones, hence the escape fraction is likely larger in 
small galaxies undergoing a burst. If this is the case, then small, 
bursting galaxies will dominate the Lyman-continuum emissivity 
even more. This strengthens our main conclusion that small, star- 
bursting galaxies can reionize the Universe by 2 ~ 10. With this in 
mind, the value of / csc = 0.2 that we used throughout this paper 

may even be conservative. 

As shown in the companion paper bv lLacev et aljj2010bh , the 
BAUGH05 GALFORM model reproduces the observed rest-frame 
1500 A luminosity function of high redshift galaxies well over the 
whole currently observed range 2 = 3 — 10. A crucial ingredient in 
this model is the boost in luminosity of galaxies as they undergo 
a minor or major merger, when the stellar initial mass function 
becomes top-heavy. This top-heavy IMF in bursts was originally 
introduced in order to fit the counts of sub-mm galaxies at much 
lower 2 ~ 2, and is necessary also to reproduce the observed high 
metallicity of gas in 2 ~ clusters of galaxies. A consequence is 
that bursts generate the majority of Lyman-continuum photons. The 
model predicts that starbursting galaxies with continuum UV mag- 
nitude Misoo.AB ~ —16, in halos of mass ~ 10 9 /i _1 M , domi- 
nate the total emissivity at 2 ~ 10 (Fig. [5}- The predict ed properties 
of these galaxies have been analysed in more detail in lLacev et al.l 
feOlObl) . Those authors show that these galaxies have stellar masses 
of M* ~ 2 x 10 5 /i _1 M n , circular velocities V c ~ 40 km s -1 , 



star formation rates M* ~ 0.06 h 1 Mq yr 1 , are gas domi- 



/Mb 



1, and have gas and stellar metallicities 



nated, M g 

of ~ 4 x 10~ a and ~ 3 x 10~ a , respectively. Assuming that on 
average approximately 2 ionizing photons are required per hydro- 
gen atom to reionize the Universe, a mean escape fraction of 20 per 
cent is sufficient to reionize the Universe by 2 = 10. 
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Figure 12. Schechter function fits (coloured lines) to the far-UV LF in the 
BAUGH05 model (heavy dots) at z = 6 over different absolute magnitude 
ranges, extrapolated to fainter luminosities (dashed lines). Starbursts (dotted 
line) introduce a deviation in the shape of the LF from a Schechter function. 
Due to this feature, varying the minimum absolute magnitude employed in 
the Schechter fit results in very different estimates of the faint-end slope 
parameter, a* (black, green and red lines correspond to minimum values 
of Afi5oo,AB of -10. -16 and -18 , respectively). Extrapolation of the fits 
to fainter values can then lead to inaccurate estimates of luminosity density 
(see Fig. ll3h which in turn results in wrongly estimated LC emissivities. 
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APPENDIX: INFERRING IONIZING EMISSIVITY FROM 
SCHECHTER FITS TO THE LF 

In observational studies, the contribution of galaxies below the de- 
tection threshold to the total ionizing emissivity is usually esti- 
mated by fitting the observed LF with a Schechter function (e.g. 
iBouwens et"aT]|2007l ; iBunker et al]|2009h . At first glance, the far- 
UV LFs predicted by the BAUGH05 model and shown in this paper 
are indeed reasonably well represented by Schechter functions, as 
they have a power-law shape at low luminosity, oc L a , and an expo- 
nential drop-off at the high luminosity end, oc exp(— L/L*). How- 
ever, the LFs predicted by GALFORM are not in detail described 
well by Schechter functions. In particular, in the BAUGH05 model, 



Figure 13. Comparison of the 1500A luminosity densities, j„, obtained by 
integrating the Schechter function fits to the BAUGH05 model LF at z = 6, 
with the values obtained from integrating over the true model LF at the same 
redshift. The colour shading indicates the ratio of the luminosity density 
from the Schechter fit to the actual value in the model. In the Schechter 
fits, M* 5 qq AB is held fixed at its best fitting value for the full luminosity 
range (black line in Fig. 1121 , while a and <f>* are the best fitting values for 
each choice of the minimum Afi5oo,AB- Values of (a, <j>*) along the dotted 
red line reproduce the actual luminosity density of the BAUGH05 model. 
Crosses mark the parameters of the three fits shown in Fig. 1121 plotted in 
the same colors. When the whole luminosity range is used for the fit (black 
cross), the model luminosity density is reasonably well reproduced by the 
integral over the Schechter function fit. On the other hand, fits performed on 
a more limited luminosity range, as in Fig. 1121 lead to significant errors in 
the luminosity density estimate. The yellow dot shows the results obtained 
with the best fit parameters from Bouwens et al. (2007) at this redshift. 



starbursts introduce a featu re (a 'bump') at ~ 2 magnitudes be- 
low at high redshifts (see lLacev et al .1 d2010bl) for more details). 
Due to this departure from the Schechter shape, the result of fitting 
a Schechter function to a GALFORM LF is strongly dependent on 
the luminosity range over which the fit is done (Fig.ll2t, Assuming 
that BAUGH05 is the 'correct' model of the high-z galaxy popula- 
tion, the observational detection limits will then strongly affect the 
predicted total ionizing emissivity, which relies on extrapolating 
the contribution of the currently unobserved low luminosity galax- 
ies based on the faint-end slope a of the Schechter fit. We want to 
investigate how much such extrapolations are likely to be in error. 

The estimated LC emissivity depends on more than just the 
LF shape, with the choice of IMF and dust extinction being cru- 
cial yet only weakly constrained by current observations. To focus 
only on the uncertainty from the assumed LF shape, in Fig.ll3lwe 
show the dependence of the 1500A luminosity density, j v , on the 
Schechter fit parameters. All values of j v were obtained by integrat- 
ing the LFs over the magnitude range —22 < Afisoo.AB < —10. 
In this figure, we vary only the normalization, <f>* , and the faint-end 
power law slope, a, and keep the characteristic absolute magnitude, 
Mi5oo,ab> fixed, because the fits shown in Fig.[T2]clearly have very 
similar M*5 o ab values. 

With this procedure, a Schechter fit over the whole luminos- 
ity range (down to Misoo.ab = —10, black line in Fig. 112b of the 
BAUGH05 LF provides a good estimate of the real luminosity den- 
sity in the model (black cross in Fig. Q~3j the luminosity density 
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from the fit is ~ 20 per cent lower than the original model). If in- 
stead the Schechter function is fit only to the brighter part of the LF 
(green and red crosses, corresponding to Misoo.AB.min of -16 and 
-18, respectively), the faint-end slope of the model is strongly over- 
estimated. As a result, the luminosity density is also overestimated 
in these cases, by factors ~ 2 and 30 for Misoo.AB.min = — 16 
and -18 respectively. We note that the BAUGH05 model predicts a 
total 1500A luminosity density at z = 6 a few times larger than 
estimates based on integrating the obse rved LF only over the cur- 
rently observed luminosity range (e.g. iBouwens et al]|2007l) . but 
this difference shrinks if the observed Schechter fits are extrapo- 
lated to lower luminosities (e.g. the yellow circle in Fig.|13|shows 
the Schechter fit found by IBouwens et all which implies a lumi- 
nosity density only 2 times lower than found in the model). 

Some authors have concluded from integrating over the ob- 
served far-UV LFs at z > 7 that galaxies alone do not emit enough 
ionizing photons to keep t he universe ionized at these redshifts 
(see e.g. lBunker et alj2009l) . but such conclusions seem premature, 
given that they do not allow for galaxies fainter than the current de- 
tection threshold or dust extinction or a different IMF slope. 

This exercise aims to point out the danger of using Schechter 
function fits to the observational data to estimate ionizing emissiv- 
ity produced by high-z galaxies. The deviations of the LF from the 
Schechter shape only add more uncertainty to the procedure which 
already hinges on a number of unknowns, e.g. the choice of the 
IMF and the dust extinction. This becomes even more important at 
higher redshifts, where the LF is even more poorly constrained by 
current observational data. 



